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. Abstract 

Iz; ■ 

' For stochastic differential equations (SDEs) with a superlinearly growing and globally 

one-sided Lipschitz continuous drift coefficient, the classical explicit Euler scheme fails 
to converge strongly to the exact solution. Recently, an explicit strongly convergent nu- 
merical scheme, called the tamed Euler method, is proposed in [Hutzenthaler, Jentzen, & 
Kloeden (2010); Strong convergence of an explicit numerical method for SDEs with non- 
'nI" , globally Lipschitz continuous coefficients, arXiv:1010.3756vl] for such SDEs. Motivated 

I by their work, we here introduce a tamed version of the Milstein scheme for SDEs with 

\^ • commutative noise. The proposed method is also explicit and easily implementable, but 

. achieves higher strong convergence order than the tamed Euler method does. In recover- 

. I ing the strong convergence order one of the new method, new difficulties arise and kind 

' of a bootstrap argument is developed to overcome them. Finally, an illustrative exam- 

ple confirms the computational efficiency of the tamed Milstein method compared to the 
tamed Euler method. 
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1 Introduction 

We consider numerical integration of stochastic differential equations (SDEs) in the Ito's sense 

dXt = fiiXt)dt + aiXt)dWt, Xo = e, te [0, T]. (1.1) 
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Here jj, : — )■ W^, a = (ui, (T2, am) '■ R-'^' — > W^^"^. We assume that Wt is an m-dimensional 
Wiener process defined on the complete probabihty space (fi, J-", P) with an increasing filtration 
{J-'t}t>o satisfying the usual conditions. And the initial data ^ is independent of the Wiener 
process. (11. ip can be interpreted mathematically as a stochastic integral equation 

ft pt 

Xt = Xo+ ^{X,)ds + y^ a,{X,)dWl t G [0,T],P-a.s., (1.2) 
^0 Jo 

where o"i(x) = {ai^i{x), ...,ad,i{x))'^ for x G G {1,2, ...,m} and the second integral is the 
Ito integral. 

This article is concerned with the strong approximation problem (see, e.g.. Section 9.3 in 
Kloeden and Platen [12]) of the SDEs (II. 2p . More precisely, on a uniform mesh with stepsize h = 
^ defined by : = to < ti < t2 < ■ ■ ■ < tp^ = T , N E N, we want to compute a numerical 
approximation Yn : M'^, n G {0, 1, N} with Yq = ^ such that 

{E\\Xt - Y^fy < ^ (1-3) 

for a given precision e > with the least possible computational effort. The strong conver- 
gence problem becomes very important because efficient Multi-Level Monte Carlo (MLMC) 
simulations rely on the strong convergence properties [5]. 

The simplest and most obvious idea to solve the strong approximation problem (II. 3p is to 
apply the explicit Euler scheme |T5] 



= K„ + hfi{Y^) + a(F„)AW„, Yo = ^,n = 0, 1, iV - 1, (1.4) 

where AWn = W^t„+i — Wt^. The method is strongly convergent with order one half if the 
coefficients fi,a satisfy the global Lipschitz condition (see, for instance, [I2]). Unfortunately, it 
has recently been shown in [7] that the explicit Euler scheme fails to provide strong convergent 
solution to the SDEs with super-linearly growing drift coefficient. It is well-known that the 
backward Euler method can promise strong convergence in this situation, see e.g.,|l]. But the 
backward Euler method is an implicit method, which requires additional computational effort 
to solve an implicit system. Recently in [9], the authors proposed an explicit method, called 
tamed Euler method, for (II. 2p 

Yr,+,=Y^ + hfi{Y^) + aiY^)AWn, with ;i(F„) = —M^^^. (1.5) 

Here /i is a modification of fi. This tamed Euler scheme is proved to converge strongly with the 
standard convergence order 0.5 to the exact solution of (11.20 if the drift coefficient function is 
globally one-sided Lipschitz continuous and has an at most polynomially growing derivative. 

On the one hand, the explicit Milstein scheme is another numerical scheme for SDEs that 
achieves a strong order of convergence higher than that of the explicit Euler scheme (II. 4p . In 
fact the explicit Milstein scheme has strong convergence order of one if the coefficient functions 
of stochastic Taylor expansions satisfy both the global Lipschitz condition and the linear growth 
condition(see [I2]). The explicit Milstein method [121 [16] applied to (II. ip reads 

m 

=F„ + M>^n) + a(F„) Aiy„ + L^'^hiyn)!-::-:^', (i.6) 

ii>i2=i 
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where 

^" = E-^..w ^nir= / dwiidw^^. (1.7) 

1^ = 1 in 'Jin 

Since the exphcit Milstein scheme and the exphcit Euler scheme coincide when apphed to the 
SDEs with additive noise, we can deduce from the results in [7] that the exphcit Milstein 
scheme generally does not converge in the mean-square sense to the exact solution solution of 
the SDEs with super-linearly growing drift coefficient. Accordingly, we follow the idea from [9] 
and replace [i{Yn) in (II. 6p with jl{Yn) to derive a tamed Milstein method 



m 



Yn+^=Yn + hjl{Yn)+a{Yn)AWn+ aj,iYn)r-^^^\ (1.8) 

Jlj2 = l 

which we expect to be strongly convergent with order one in the non-globally Lipschitz case. 

On the other hand, although Milstein-type schemes may achieve a strong convergence order 
higher than that of Euler type schemes, additional computational effort is required to approxi- 
mate the iterated Ito integrals for every time step [13]. This will enable the Milstein-type 
schemes to lose their advantage over the Euler-type schemes in computational efficiency. In this 
article we restrict our attention to SDEs with commutative noise, in which case the Milstein 
scheme can be easily implemented without simulating the iterated Ito integrals. In this sit- 
uation, Milstein-type method is much more computational efficient than Euler-type method. 
More precisely, let the diffusion matrix a fulfill the so-called commutativity condition: 

L^^(^k,j2 = L^^(^k,n, ji,j2 = l,---,m,k = l,...,d. (1.9) 

In many applications the considered SDE systems possess commutative noise (see [T2]). 

Thanks to the property + 4"^^' = ^W^^AWl^Ji ^ ja, in this case the tamed 

Milstein method ( II. 8p takes a simple form as 

^ m 

y„+i =y„ + hfiiYn) + aiY^)AWn + 2 5Z ^''^^^l^^n) {AW^^AW^^ - 6j,j,h) , (l.io) 

ii,i2=i 

where Sj^^j^ = 1 for ji = j2 and Sj^^j^ = for ji 7^ ja, /i is the modification of /i as defined in 

The main result of this article shows that the tamed Milstein scheme ( ll.lOp converges 
strongly with the standard convergence order one to the exact solution of SDEs with com- 
mutative noise if the drift coefficient fi is globally one-sided Lipschitz continuous and has at 
most polynomially growing first and second derivatives. The diffusion coefficient a and the 
coefficient function of stochastic Taylor expansion L^^aj^, ji, j2 G {1, ...,m} are assumed to be 
globally Lipschitz continuous. It is worthwhile to mention that a similar approach as used in [9] 
is evoked to obtain uniform boundedness of p-th moments of numerical solutions produced by 
the tamed Milstein method. We also introduce similar stochastic processes Dn that dominate 
the tamed Milstein approximation on appropriate subevents fin (see Section 2 for more details). 
With bounded p-th moments at hand, our main effort is to show for the time continuous tamed 
Milstein methods there exists a family of real numbers Cp^T G [1, 00) for p e [1, 00) such that 

^^^<Cp,T-h, he (0,1]- (1-11) 





sup 


Xt-Yt 


p- 


) 




^te[o,T] 
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The key difficulty is that relative to previous analysis [9] a sharper estimate of the term J (see 
f l3.35p ) must be obtained to get the strong convergence order one. To overcome this difficulty, a 
certain kind of bootstrap argument is exploited (see the estimate of J for more details). To the 
best of our knowledge, this is the very first paper to successfully recover the strong convergence 
order one for the Milstein-type method under non-globally Lipschitz condition. 

The rest of this paper is arranged as follows. In the next section, uniform boundedness of p-th 
moments are obtained. And then the strong convergence order of the tamed Milstein method is 
established in Section 3. Finally, an illustrative example conffims the strong convergence order 
of one and the computational efficiency of this scheme compared to the tamed Euler scheme. 



2 Uniform boundedness of ^^-th moments 

Throughout this article, G N is the step number of the uniform mesh defined in the previous 
section. Moreover, we use the notation ||x|| := (|xip + ... + 1x^^)2, {x,y) := Xiyi + ... + x^yk 
for all X = {xi,X2,...,Xk),y = (l/i, I/2, eR'',k E N, and \\A\\ := sup^gjji ||^||<i \\Ax\\ for all 

A G M'^^', /c, / G M. Furthermore, we make the following assumptions. 

Assumption 2.1 Let fi{x) and ai{x),i = 1, ...,m be continuously differentiahle and there exist 
positive constants K >1 and c > 1, such that Vx, y 

{x-y,n{x) - niy)) < KHx-yf, (2.1) 

Mx)-a{y)\\ < K\\x-yl (2.2) 

\\U'a,,{x)-U'ajM\\ < Ji,j2 G {l,2,...,m}, (2.3) 

y{x)\\ < K{l + \\xf). (2.4) 

Note that the globally one-sided Lipschitz condition (12. ip on the drift /i and the globally 
Lipschitz condition (12.21) on the diffusion a have been widely used in the literatures, e.g., 

[HEIEIITIIHIE]. 

To prove uniform boundedness of p-th moments of the numerical solution, we follow the 
ideas in [9] to introduce the appropriate subevents fi„ and dominating stochastic processes 

fi„:=|a;Gfi| sup Dk{uj) < ^ sup ||AW^fc|| < l|, (2.5) 

0<fc<n-l 0<fc<n-l ^ 

n-1 

Z}„:=(A+||e||)exp(A+ sup V [A|| AW^^f + a,] ) , (2.6) 



0<u<n 
— — k=u 



where 
A = 
and 



4TK + 2Ty{<d)\\+2K + 2\\(y{Q)\\+m^{T +1)(k + max \\U'ajM\\ 



Oin -•iii/„ii-'ifxi|-^||7 ' ^ll|-'n||i-l-J-\||-I^||7 / J OWV 

JiJ2=l 

(2.7) 

The following lemmas are needed in order to prove uniform boundedness of p-th moments. 
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Lemma 2.2 Let Yn,Dn and f2„ be given by U . 1 0\) . I[27^) and Ii2. 5\) . respectively. Then 

<Dn, for all n = 0,1,..., N. (2.8) 

Proof. First of all, note that \\AWn\\ < 1 on for all < n < - 1 and N e N. The 

globally Lipschitz continuity of cr and U^aj^, ji, j2 £ {1,2, ...,m}, and the polynomial growth 
bound on /i' imply that, on fln+i H {u; G < 1} for all < n < — 1, 



n+1 1 



ii,i2=i 

<1 + - /i(0)|| + h\\fiiO)\\ + ||(T(r„) - (t(0)||||AW„|| + ||(t(0)||||AW„|| 

h,h=i (2.9) 

<i + hK{i + llFnDllFnll + /i||Mo)|| + i^||y;||||Aw^„|| + ||a(o)||||Aiy„|| 

+ 2 E (i^rn|| + ||L^V,,(0)||)|Aiy;j^Aiy^-^-5,,,,/i| 

<1 + 2TK + T||/.(0)|| + + ||a(0)|| + -(m + T) [k + max ||L^^a,,(0)|| 
<A. 

2 l2 

Moreover, the Cauchy-Schwarz inequality and the inequality a ■ h < ^ + ^ ioi all a,h E M. give 
that 

||F„+if = \\Yn + /i/i(r„) + a(y„)Aiy„ + M„||' 

2 I !,2|| r-./-^/- M|2 , ll_/AZ\ATJ/ l|2 , Ilji/T l|2 



/i^llMi;)lr + lk(>;^)Aiy„ll^ + l|M„ 
+ 2/i(r„, /i(r„)) + 2(r„, a{Y^)AWn) + 2(r„, m„) 

+ 2h{~^{Yr,), a{Yn)AWn) + 2/i(/i(yO, M„) + 2{a{Yn)AWn, M„) (2.10) 
<||r„f + ?,h^MYr,)f + 3||a(yOin|AIV„|P + 3||MJ|2 
2/i 

/i(K)) + 2(K, a(FO AVTn) + 2(F„, M„) 



l + hMY^)\\ 
on f2 and < < — 1. Here we denote 

Mn = ^Y. L^'^niyn){AW^^AW.i- - 5,,,,/i). (2.11) 

ii J2=i 

Additionally, the global Lipschitz continuity of a,U^aj2 implies that for > 1 

\\a{x)r < ilWix) - am + ||a(0)||)^ < {K\\x\\ + ||a(0)||)^ < {K + \\am)W, (2-12) 

and 

||L^V,,(a:)f < (||L^V,,(x) - L^V,,(0)|| + ||L^V,,(0)||)2 

<(Ji||x|| + ||L^V,,(0)||)2 (2.13) 
<{K+W^a,M\\)W, 
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and the globally one-sided Lipschitz continuity of /i gives that 

= -/i(0)) + (x,/i(0)) < A'llxf + ||x|| ■ ||/i(0)|| < (A^+ ||/i(0)||)||xf . (2.14) 

Furthermore, the polynomial growth bound on /i' implies that 

||/i(x)f < (Mx) - m\\ + IIMO)||)' < mi + ||xr)||a:|| + M0)\\)" 



\ic+l) 



wmwy 



|2(c+l) 



< (2/^11 

<(2A^+||/i(0)||f ||x 
<iV(2K+||M0)||f 

on 1 < < iV^. Combining f l232|) - fl2:T5|) . we get from i^AO\i that 



\\Yn+if <r„f + 3hT {2K + \\^^mf ^nf + 3(A^ + ||a(0) ||)lr„f || AP^^f 
+ 3|| A4f + 2/i(A' + ||/i(0)||)||F„f + 2(F„, a(K)AW^„) + 2(r„, M„) 



(2.15) 



(2.16) 



on {oj e VL\l < \\Yn{uj)\\ < N^}. Since u G Qn+i implies ||APV„|| < 1, on cu G Qn+i n{a; G |1 < 
||>"n(c^)|| < N^c}^ we derive from (12TT]) and (|27[3|) that 



Jl J2 = l 
2 »^ 



J1J2=1 

,3 , ,2 



m 



< — {K+ max \\L^'a 



n 



|2 

-nil ; 



2 V i<ii,i2<"i 

where the fact that lAl^^^p < ||Aiy„f < 1 was used. Inserting fl2Tfj) into fl2A6l) shows that 

||^n+lir 

<\\ynf [l + 3^(2i^ + Mm? + 3(A' + ||a(0)||)l Aiy„|r + 2^(1^ + ||/i(0)||) 



+ max ||L^V,,(0)||) (||AW^„f + /i^ 



< 



2 /3 



V2 



(2Ti^ + T||M0)||)^ + T(A'+||M0)||) 



2 + 2(F„, a(FO Aiy„) + 2(F„, 

V+ max \\U^a,M\\) 



4 



+ 2(|(A^+||a(0)||)^ + ^ 



(k+ max ||L^V,,(0)||)'^ ||AW^„||2 + 2a„ 

V l<ii,i2<m / / 



< 



exp 



2A 

iV 



+ 2A||Aiy„||2 + 2a„ 



(2.18) 



on f2„+in{a; G |1 < ||y„(aj)|| < A^2c}. Now combining (12. 9 p and (12.18^ . and using mathematical 
induction method as used in the proof of Lemma 2.1 in [9] finish the proof. □ 
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Lemma 2.3 For all p > 1 



sup E 

N>A\pT 



exp (p\^\\/^Wk\ 



< oo. 



(2.19) 



k=0 



Proof. This result is identical to Lemma 2.3 in [9] with only different A. □ 
Lemma 2.4 Let an be given by \2. 7\). Then for all p > 1 



sup sup E 

z£{-l,l} 



n-1 

sup exp I pz 2, oik 

0<n<N V 



< OO. 



(2.20) 



Proof. We set 



1 1 / o-{Yk) 



2 I ^ h V — ^ 



iij2=i 



Then we have = + a^. Hence Holder's inequality shows that 



n-l 



sup exp [ z ftfc 



n-l 



fc=0 



< 



sup exp [z V 

0<n<A' ^ — 



n-l 



fc=0 



L2p(f7;: 



sup exp [z^al 

0<n<N ^ — 



fc=0 



Note that Lemma 2.4 in [9] has proved that 



sup 

2e{-i,i} 



sup exp [ z 

0<n<N 



n-l 
k=0 



L2p(f7;] 



< OO. 



L2p(n;M 
(2.21) 



(2.22) 



Consequently it remains to prove the boundedness of the second term on the right-hand side 
of fl2.2ip . One can easily verify that the discrete stochastic process z Yl^Zo ^t^'"^ ^ {0, 1, N} 
is an {J^tn : < n < A^}-martingale for every z G { — 1, 1}. Since the exponential function is 
convex, the discrete stochastic process exp (2; ^^Zq a^) ,n G {0,1,. ..,A^} is a positive {J^tn '■ 
< n < A^}-submartingale for every z G {—1,1}. Therefore, the Doob's maximal inequality 
gives that 



n-l 

sup exp I z ot\ 

0<n<N V 



< 



2p 



2p-l 



N-l 

exp ( 2: ^ al 

k=0 



L2P{Q;I 



(2.23) 



The Cauchy-Schwarz inequality, f l2.13p and Lemma 12.31 give that 



7V-1 

exp ( 2 ^ al 

k=0 



L2p(C;] 



< 



k=0 



ji,j2 = l 



2im 



N-1 



< 



< 



f X K 1 

exp ( > ( 1 — max \\U^(y 

fe=o ji,i2=i 



5^ (jAW^^^Aiyf 

L,i2 = l 

V — ' I I l<Jij2<m V 
fc=0 

N-\ 

exp(A5^||Aiy,| 



^il,i2^ 



L2P{n;] 



L2P(0;1 



< OO. 



L2P(Q;] 



This together with fl2.23p completes the proof. □ 

Lemma 2.5 Let Dn he given by Ii2.6\) . Then for all N > 8XpT and p > 1 



sup E 



sup \Dn\^ 

0<n<N 



< OO. 



(2.24) 



Proof. Note that Dn here takes the same form as in with only different With 
Lemma [2.31 and Lemma [2.41 at hand, one can follow the proof of Lemma 2.5 in [9] to derive the 
desired result. □ 



Lemma 2.6 Let Qn be given by Ii2.5\) with n = N. Then for all p > 1 

sup {N^FlinNT]) < OO. 

TVgN 



(2.25) 



Proof. The proof is identical to the proof of Lemma 2.6 in [9]. □ 

Before establishing the main result of this section, we also need two Burkholder-Davis-Gundy 
type inequalities. 

Lemma 2.7 Let k E N and let Z : [0,T] x f2 — )■ be a predictable stochastic process 

satisfying ^{Jq WZgW^ds < oo) = 1. Then for all t G [0, T] and all p > 2 



sup 

se[o,t] 



ZudW^ 



< p 



t "I 



i=l 

Here the vectors Ci = (1,0, ...,0)^ G W''\ eg = (0, 1,...,0)^ G M™, 
are orthogonal basis of vector space M"^ . 



(2.26) 
= (0,0,...,!)^ G M"^ 



Proof. Combining Doob's maximal inequality and Lemma 7.7 of Da Prato, G., and Zabczyk 
[2] gives the desired assertion. □ 

The following is a discrete version of the Burkholder-Davis-Gundy type inequality (12.261) . 

Lemma 2.8 Let k E N and let Zi : Vt ^ M'^^'",/ G {0, 1, ...,N — 1} be a family of mappings 
such that Zi is J-'it/B{R^^"^) -measurable. Then for all < n < N and p > 2 



sup WJ^Zi^Wi 



n—l m 



0<j<n 



1=0 



LP{n;] 



1=0 i=l 



T 

«llLP(n;R'=) jY 



(2.27) 
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Theorem 2.9 Let Yn be given by U.10\) . Then for all p G [1, oo) 



sup 



sup E||F„r 

0<n<Ar 



< OO. 



(2.28) 



Proof. From ( ll.lOp we have 

ll^n||LP(f7;K<^) ^ II 'Cll LP(n;lR'') 



n-1 



n-1 



n-1 



k=0 



LP{n;] 



fc=0 



LP(0;] 



k=0 



(2.29) 



LP(Q;I 



where the notation comes from 02.111) . Using (12.21) . the triangle inequahty and the Burkholder- 
Davis-Gundy type inequahty in Lemma [2.81 we have 



n-1 



n-1 



n-1 



fc=0 



LP(Q;1 



< V [a{Yk) - a(0)] AWk + V a{0)AW, 



fc=0 
n— 1 m 



fc=0 



lLP(n;IE 



1/2 TlT -^/^ 

^^'l E E - ^^(o)iil(n;M^) ^) + E 11^^(0)11' ) 

fc=0 j=l i=l 



1/2 

<p(mi^2^^||n,||J,(^^j,,)) +pv^||a(0)||. 



A:=0 



(2.30) 



For the fourth term on the right-hand side of (I2.29p . the estimate in the second inequahty of 
(I2.13P and the independence of Yk and AWk imply that 



71, — 'L Tfh 



/c=0 

n—1 m 



k=0 ji,i2=l 



LP(Q;I 



(2.3i; 



lJ2^ILp(n;] 



k=0 jl J2=l 

Using Burkholder-Davis-Gundy type inequality (I2.27P and mutual independence of AWk gives 

Jl,i2 = l jlJ2 = l 

=^ E l|AW-f|U.(n;M)-||Aiyf|U.(n;M) + -^||AIV^||i2.(o^^^ ^2.32) 



ii J2=i 
ii7^i2 



where c 



rt-1 

n 1 

k 



p2 + 2mp2 + ^. Inserting fl232|) into f l23T]) we obtain that 



E^^ 



n-1 



fc=0 



^ < /"^Cp^r^^V" II ifc II Lp(n;i8d) +Tcp,m sup \\L^^aj^{0)\\. 

LP{n;R'i) ^ l<ii,j2<m 



(2.33) 
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Combining f l2:29|) .f l230|) and f l233|) . and using \\fi(Yk)\\ < 1 give 

1/2 

ll^n||Lp(f7;Kd) < 11^ llLp(f7;Md) + N + p(mK^h ^ || ^fc' || L(n;Md) ) + pVfm]\a(0)\\ 



k=0 

n-1 

+ Kcp^mh^\\Yk\\ LP (n;R'i) +Tcp^m sup \\L^^ aj^{0)\\. 

1<J1 J2<m 



(2.34) 



fc=0 

Thus taking square of both sides shows that 



<^{U\\Lv{n-,Rd) + N + pVTm\\a{0)\\+Tcp^rn sup \\U'aj^{0)\\ 

1<JU'2<™ 

n— 1 rt— 1 

fc=0 fc=0 



(2.35) 



In the next step Gronwall's lemma gives that 

sup ||F„||ip(f^,K'*) < V3 exp(Ci) (||{||LP(n;iRd) +N + C2) , (2.36) 

^ AT ^ ^ 



0<n<N 



where Ci = l{3mKY + 3K^cl^T)T, C2 = pVTni\\a{0)\\ + Tcp,mSup^<j,j,<m\\L''aj,{0)\\- 
Due to the on the right-hand side of 02.361) . f l2.36p does not complete the prove. However, 
exploiting f l2.36p in an appropriate bootstrap argument will enable us to establish f l2.28p . First, 
Holder's inequality. Lemma [2.61 and the estimate ( I2.36P show that 

sup sup \\l{n,,rYn\\ LP in-M'i) 

7VGNO<n<7V 

<SUp sup (||l(n„)'=||L2p(n;K<*)||^n||L2p(n;Rrf)) 
7VGNO<n<7V 



<( sup (iV- ||l(f^^)c||^2p(t^.Kd)) )( sup sup {N ■ \\Yn\\L2p(n;R''))) (2.37) 
^A^GN ^ ^ NmO<n<N ^ 

<(supiV2*^-P[(l];v)l)'^(sup sup (iV-^-||F„|U2.(^;i,.)) 



■iVGN ' ^ NGNO<n<N 

<oo. 

In addition. Lemma 12.21 and Lemma 12.51 imply that 

sup sup ||ln„^n||LP{C;IRd) < sup SUp ||D„||ip(f^.Kd) < CX). (2.38) 

NG^O<n<N NeNO<n<N 

Combining f l2.37p and fl2.38p finally completes the proof. □ 



3 Strong convergence order of the tamed Milstein method 

In order to recover the strong convergence order one for the tamed Milstein method, we ad- 
ditionally need the following assumptions. Throughout this section Cp^T is a generic constant 
that might vary from one place to another and depends on /i, a, the initial data ^, and the 
interval of integration [0,T], but is independent of the discretisation parameter. 
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Assumption 3.1 Assume that fi{x) andai{x) are two times continuous differentiahle and there 
exist positive constants K,q > 1, such that Vx G W^,i = 1, ...,d 



l/^"(a;)||L(2)(Rd;E<*) < K{l + \\xp) 

Wii^)\\L(^)iRd;R'i) < ^• 



(3.1) 
(3.2) 



Here, for a two times continuous differentiahle function / : M*^ — > R*^ we use the notation 
\\f"ix)\\L(2){R'';R'') = ^^Ph,,h2eR'',\\h,\\<i,\\h2\\<i ll/"(a;)(^i,^2)||. The bilinear operator f"{x) : M'^ x 
R'^ M"' is defined by 03.81) . In the following convergence analysis, we prefer to write the tamed 
Milstein method in the form of fll.81) . rather than fll.lOl] . We now introduce appropriate time 
continuous interpolations of the time discrete numerical approximations. More accurately, we 
define the time continuous approximation Ys such that for s G [tn,tn+i) 



iij2=i 



/g III ng Ttl fig TTi 

li{Y^)du + Y. <yn)dwi + Y. Xl^v,(F„)Aiy^ciiy: 



(3.3) 



where we use the notation 



It is evident that Yt^ = Yn,n = 0,1,.. .,N, that is, Yt coincides with the discrete solutions at 
the grid-points. We can rewrite Yt as an integral form in the whole interval [0, T] 



i=l 



%=Yo+ I jl{Y^Ms + J2 I cr.(r„j + 5^LV.(y„jAiy; 







dWl 



(3.4) 



where Ug is the greatest integer number such that tn, < s. Combining (11.21) and (13. 4p gives 

m „4 m 



i=l 



Xt-Yt= I [/i(X,) - /i(r„j] ds + J2 I [^^iXs) - a.(y„j - J2 L'cr,{Y^J^Wi 



dWl 
(3.5) 



In what follows, we also use deterministic Taylor formula frequently. If a function / : M*^ — 
is twice differentiahle, the following Taylor's formula holds 



/(n) - /(F„J = /'(F„J(n - F„J + i?i(/), 
where Ri{f) is the remainder term 

Riif) = f\l - r)f"{Yn, + r{Z - YnMYs " Y,,^,Y, - F„Jdr. 
Jo 



(3.6) 



(3.7) 
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Here for arbitrary a, hi,h2 G M.'^ the derivatives have the following expression 



(3i 



Replacing — in (13.61) with (13. Sp and rearranging lead to 
where 



(3.9) 



m 

Riif) = nynj{is-t^Mynj+ Yl L^'^n(Ynji^:-:)+Ri{f). (3.10) 



ii,i2=i 

By the definitions (11.71) and (13. 8p . it can be readily checked that 
Therefore replacing / in (13. 9p by (Xj and taking (13.111) into account show that 

m 



(3.11) 



(3.12) 



Theorem 3.2 Let conditions in Assumptions \2. 1\ and \3.1\ and he fulfilled. Then there 

exists a family of real numbers Cp T > 1, P > 1 such that 



i/p 



E 



sup \\Xt-n 



<Cp^T-K /iG(0, 1]. 



(3.13) 



The following five lemmas are used in the proof of Theorem 13. 2[ 



Lemma 3.3 Let conditions in Assumptions \2. 1\ be fulfilled. Then for all p > 1,1 < ji, j2 < 
the following estimates hold 



sup sup 

Nef>SO<n<N 



E|iMF^)rV^ii/^'(^-)ii'V^ii^(^")ii'V^ii^''^^-^(^") 



< oo. 



(3.14) 



Proof. It immediately follows from Theorem 12.91 by considering (12. 4p and (12. 121) - (12. 151) . □ 
Lemma 3.4 Let conditions in Theorem \3.2\ be fulfilled. Then for all p > 1 



sup 

te[o,T] 



\Xt\\LP{n;R'') V \\K^t)\\LP{n;V) \/ IWi^i 



tj|Up(Q;M'*x'") 



t\\LP{n;] 



< OO. (3.15) 



Proof Conditions (I2.ip - (l2.2p ensure that the exact solution Xt satisfies sup^g^.T] ^^11^* P < ^ 
(See, e.g., [TU Theorem 4.1]). Polynomial growth condition on /x as (I2.15P and linear growth 
condition on a give the second and the third estimates. Taking the definition (13. 3p and Lemma 
13.31 into account, we can easily obtain the last estimate. □ 
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Lemma 3.5 Let conditions in Theorem \3.2\ be fulfilled. Then there exists a family of constants 
Cp^T > 1 such that for p > 1 



sup \\Yt-Yn,\\LP{n;Rd)\/ SUp \\fi{Xt) - fi{Xt^J\\LP{n;] 
tG[0,T] t£[0,T] 



<Cp,Th^. (3.16) 
Proof. Let nt be the greatest integer number such that t„j < t. From (13. 3p we have 

(3.17) 



%-Y^, = {t-tMYnJ + a{YnMWt-WtJ+ J2 L''^j,iYnJI^:-l 

ii J2=i 

Following the same line as estimating (I2.29p . one can readily derive the first estimate. For the 
second estimate, we use (12. 4p and Holder's inequality to obtain 



IIM^*) - ^(x,„j 



<K\\(i + \\Xtr + \\Xt\n-\\Xt-x. 



<Ai|(i + iix,r+iix r) 



■\\Xt-Xt 

1 1 ^ "-'It I 



L2p(Q;] 



where 



^^{X,)ds+ / <y{Xs)dw,. 



(3.18) 



(3.19) 



Taking Lemma 13.41 into account and using the same argument as in the previous section, one 



can obtain \\Xt ~ Xt„^ Wl'^p^q-,] 
flXTSj) and Lemma |331 □ 



< Cp^rh^ , and then complete the proof easily by considering 



Lemma 3.6 Let p > I and let conditions in Theorem \3.2\ be fulfilled. Then for i = 1,2, m 

II -Ri(/^) II Lp(n;Md) V \\Ri{cri)\\Lp{n;M.'i) < Cp^rh. (3.20) 

Proof To estimate \\Ri{fi)\\LP(n;R'i), we need to estimate ||-Ri(/i) ||LP(n;iR'*) first. Due to Theorem 
12.91 and Lemma 13.41 we can find some suitable constant Cp^x such that 



LP(n;IE 



< / (1-0 
•1 

< 



dr 



LP{n;] 







(Yns + r{Ys - >"nJ)||L(2)(Rd.Kd) ■ ||1". - Yn 

<i^||(i + ||F„j|'' + ||F,r)-||n-i; "211 

Ill + UK, 11''+ IlKir'" 



LP(Q;1 



dr 



(3.21) 



lLP(Q;I 



lL2p(f7;l 



Y —Y 



<Cp^Th, 



where the polynomial growth condition (13.11) on fi"{x). Lemma 13.51 Holder's inequality and 
Jensen's inequality were also used. Now we return to ||-Ri(/i)||Lp(n;R'*)- Replacing / in (I3.10p 
with /i gives 



\Rl if') \\ LP (n;R'i} 



(3.22) 
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iij2=i 



LP(t7;R'*) 



l^l(/i)ll 



Lp{n]] 



^ m 



ii J2=i 



LP(n;Rrf) 



LP(n;IR'i) • 



Following the same line as estimating (I2.32p and noticing that s — < h, one can similarly 
arrive at 

rn 

- \\^wi^^wi^-K,nHLnnm<(^v,mh. (3.23) 

ji,i2=i 

Further, using Holder's inequality we derive from Lemma 13.31 that for < s < t < T, 1 < 

ji, j2 < rn 



ll/^'(^nJ/^(^nJ||LP(C;R'*) < II /^'(K J || L2p(Q.Kdxd) ■ || /i(^r!. J || L2p(n;Rd) < 



OO, 



|/i'(F„jL^-V,,(F„, 



lLP(n;I8'*) 



< llAt'('^nJ||L2p(f^.IKdxd) ■ IlL^Vj^CKriJl 



L2p(n;IE 



< OO. 



(3.24) 



Now, using the independence of and Aiy/s AW/^ and combining fl3:2TD . fl3:22|) . fl3:23|) and 
(I3.24p . one can show 

ll^i(/^)l|Lp(n;Rd) < Cp^rh. (3.25) 
In a similar way as estimating ||_Ri(/i)| 



LP(n;R<*)! 



one can derive that 



|-Rl(<7j) ||LP(n;Kd) < Cp^rh. □ 



(3.26) 



Our proof of Theorem 13.21 also needs the following Burkholder-Davis-Gundy type inequality for 
discrete-time martingale (see Theorem 3.28 in [H] and Lemma 5.1 in [8]). 

Lemma 3.7 Let Zi,...,Zn : Q ^ be T/BiM.) -measurable mappings with K\\Zn\\^ < oo for 
all n e {1, A^} and with E[Z'„+i|Z'i, Z„] = for all n 6 {1, A^}. Then 

\\Zi + ... + ^„,||LP(n;R) < Cp (||-Z'i|||p(Q.]g) + ... + ||^n||LP(Q;M)) ^ (3.27) 

for every p G [2, oo), where Cp are constants dependent of p, but independent of n. 
Proof of Theorem \3.2[ Applying Ito's formula to f l3.5p gives 

||X, - =2^ (X„ - F„,/i(X„) - ftiYnj)du 

m „g m 

+ 2J2 (x„-F„,a,(X„)-a,(K„J-^LV,(F„JAW^^)diy^ .g^g) 
i=i j=i 

m „g m 2 

+ ^ / ILiX^) - a,(F„J - J2 L'a.iY^JAWl du. 
i=iJ^ " ,=1 
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For the integrand of the first term in fl3.28p . we use (12. ip and the Cauchy-Schwarz inequahty 
to arrive at 



X„-F,,/i||/i(F„J||/i(F„ 



(3.29) 



For the integrand of the third term in 
(13:T2|) and (Q to get 



one can use an elementary inequahty, the notation 



a,iX^) - (T,(F„J - J2 L'a,{Y^J^Wl\\ <2||a,(Xj - a,(Fjf + 2||i?i(a,) 

<2i^'||X„-Kf + 2||i?i(a,)f. 
Inserting (I3.29P and (I3.30p into (I3.28P and using the notation (13.121) show 

PS PS 

WX.-YsW' ={2K + l + 2mK^) \\Xu - Y^fdu + ||M>"„Jfrfn 



(3.30) 



+ 2^/ \\RMWdu + 2 (x,-K,/i(K)-/i(F„J)rfM 
Jo Jo \ 

+ 2 ^ / (X„ - F„, (T,(X„) - a,(Fj + i?i(a,))rfiy:. 
i=i 



(3.31) 



Hence 



sup \\X,-Y,\\' <{2K + l + 2mK') / \\Xs -Y^U'ds + h' / ||/i(F„Jfrfs 

0<s<t 



+ 2V/ ||i?i(a,)fds + 2 sup / (X„-F„,/i(Fj-/i(F„J)rfM 

~7J0 0<s<iJo \ ' 

m „s 

+ 2 sup V / (X„ - F„, a,(X„) - a,(K) + R,{a,))dWi 

0<s<t~^ Jo \ I 



and thus for p > 4 



sup ||Xs - F 

0<s<t 



LP(n;M) 



sup ||X^ - Ys\ 



0<s<t 



p 



<{2K + 1 + 2mK') [ \\Xs-Z\\l^^.^^,)ds + h' [ MY^MUin;R-)ds 
Jo Jo 



sup 

0<s<t Jo 



Xu - ^^{Yu) - fJ^iYnJ )du 



+ 2J2 \\RM 

i=l -^0 

m „g 

+ 2|| sup V / (X„ - F„, a,(Xj - a,(Fj + R^{a,))dW: 



(3.32) 



L7(C;M) 
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For the last term on the right-hand side of (13.321) . we use (12. 2p . (I2.26p . the Cauchy-Schwarz 
inequahty and elementary inequalities to derive 



M))dw: 



2 sup V / (Xu- aiiXu) - (TiiYu) + R 

" 0<s<t ~Y "'0 ^ 
m 

< 2V||sup / (X„-K,a,(X„)-a,(K) + i?i(a,))ciW^: 

~^ II 0<s<f Jo 



< 



P 



E 

i=l 



ds 



i=i ^-^0 " 



< snp \\Xs -Ys\\LP(^n;R'i) ■ P^i \\(ri{Xs) - ai{Ys) + Ri{ai) 



0<s<t 



1=1 



LP(n;R<*) 
2 



ds 



< - sup - +p2^y^ / - ai{Ys) + 

+ 2p^m^K [ \\Xs-Z\\l,^^.^,^ds 
Jo 



LP{n;R'') 
2 



ds 



ds 



< \ sup - Y, 

+2p2mV / \\RM. 
i=i -^0 



(3.33) 



At the same time, replacing / in (13. 9 p by yU and using the Cauchy-Schwarz inequality and 
elementary inequalities give 



2 

=2 
<J + 2 



sup / (X„ - F^.,/x(y„) - /i(F„^))rfn 

0<s<tJo \ ' 



sup / (X„-K,/i'(r„J(a(K„JAl^„) +i2i(/i))rfM 

0<s<tJo \ I 



p 



V 



sup / (Xu -Yu,Ri{fi))du 

0<s<tJo \ I 

t 

2 J„ I / II 5 /,.M|2 



Jo Jo 



where we denote 



J = 2 



sup / (Xu-Yu,fi'{YnJa{YnJAWu)du 

0<s<tJo \ ' 



L^(C;1R) 



(3.34) 



(3.35) 
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Inserting (K33^ and flSlMj) to f l332|) yields 



sup ||X^ - Ys\ 

0<s<t 



LP{n;R) 



<2{K + l + mK'+p'm'K) \\Xs -Y,\\i,^^.^,^ds + / MYnM^.^nm^^ (3.36) 



+ 2(1 +p^m) 



)\\LP{n;R'i-)^^ + / ||-Rl(/i)||^p(-Q.igd)(is + J. 



Therefore it remains to estimate J as (13.351) . Using (13.31) . (I3.19P and (I3.12p shows 



i = l "-nil 



+ {u- t„j /i(Xi„j - (m - t„j /i(r„j + Xt^^ - Yn^. 



Thus 



J <2 



sup 



[/.(X,) - /i(Xi„J]rfr,/i'(F„Ja(r„JAiy„)rf« 



L^(C;IR) 



0<s<t Jo 



i = l '"u 



+ 2 

+ 2 



/ (E/ ^i(^.)^^W^;,/i'(>^nJa(y.JAW„\dn 
sup / (c„,„,/i'(F„,J(T(F„JAlK)f^M p 

)<s<tJo \ I L^(n;R) 

sup r(x,„„ -F„„,/i'(F„Ja(r„JAiy„)rf« 



0<s<t 

= Jl + J2 + ^3 + -''4 + ^5 



L^(C;1R) 



where Cn„ e is defined by 

Cn„ = (m - t„J /i(Xt„J - (m - t„J /i(FnJ. 

To begin with, we estabhsh the estimate 

||/i'(F„Ja(r„JAiy„||ip(f^.Kd) < C^prhh, 



(3.37) 
(3.38) 
(3.39) 



which is frequently used later. Using Holder's inequality, Lemma [3.31 and Lemma [2771 we know 
that for < A; < X - 1, tfc < s < 4+1 



||/^'('i^fc)0"(^fc)|lL2p(j^.]Kdxm) < ||/i'(^fc)||L4p(n;Rdxd)lk('^fc)||L4p(n;Rd>. 



< 00 



(3.40) 
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and 



{3A1] 



Combining f l3.40p and fl3.4ip one can readily obtain f l3.39p by Holder's inequality. Concerning 
Ji, we use Holder's inequality, fl3.16p and fl3.39p to arrive at 



Ji < 2 



t fU 



Jt, 



^^'{Yn^)aiYn^)/\Wu drdu < C^^rh^ (3.42) 



For J2, using Holder's inequality, (12. 2p . (12. 26 p . elementary inequalities and (13.390 gives 



" i=i -'w 



i = l >' l-nu 

•2 H ru ™ 



LP(n;I8d) 
2 

LP{Q;1 



LP(n;IE 



du + h 



LP{n;l 



du 



2 ft f"U '"' 

^ T $^||a,(X.)-a,(n: 



< mp^K'^ / sup — Fr"^ 

Jo 0<r<u 



LP(n;I8<*) 



(3.43) 



For J3, similarly as above we obtain that 

ft ,, r« _ 

i?l((Ti)rfWC 

LP(n;l 



/" II . _ 

J3<2 / 5^ / RM)dWl 
Jo II 

m 



<2p 



^ >/t 



dr 



f^'{Y^Ja{Y,JAW. 

LP{Q 

^^'{Y^Ja{Y^JAWu 



du 



du 



(3.44) 



LP{Q;R'') 



<Cr,Th^ 



where (I3.20p and (13.390 were also used. Now, it remains to estimate J4 and J5. We split J4 
into two terms as follows: 



Ja <2 



n.-l 



+ 2 



U,^i'{Y^MynJAW^du 



sup 

0<s<t I 



L7(C;1 



(3.45) 



:— J41 + Ji2- 

Recall that (k ^ J^t^: for = 0, 1, — 1. It can be readily verified that the discrete time 
process 



J2 ['^\Ck,f^'iYk)aiYk)AW^yu\, ne{0,l,...,iV} 
. fc=0 "^'fc ' 
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is an {Ttn : < < A^} -martingale. With the aid of Doob's maximal inequality, Lemma 13.71 
and Holder's inequality we obtain that for p > 4 



J41 < 



< 



2p I 
p-2\ 



a,/i'(n)a(n)AiyJciM 



p-2 



(5: 



fc=0 



tk 



Ck,fi'{Y,)a{Yk)AWu)du 



L7(n;R) 



1/2 



< 



< 



f fc=0 



fc=0 



a,/i'(n)a(n)Aiy„ 



1/2 



y fc=o 



1/2 



where by fl3.38p , Lemma 13.31 and Lemma 13.41 

IK*^llLP{n;M'*) - Hf^iXtJ\\LP{n;K'') + h\\Kyk)\\LP{n;K'i) < Cp^xh. 

Hence, taking fl3.39p and (13.471) into account, we derive from fl3.46p that 
For the second term J42, Holder's inequality, fl3.47p and f l3.39p give that for p > 4 

'J42 \ 2 

T 



(3.46) 



(3.47) 



(3.48) 



=E ( sup 

■ 0<s<t 



<h2-'E{ sup 

■o<s<tJtn, 



(C„.,/i'(i;ja(i;jAiy„) 

{Ck,l^'{Yk)a{Yk)AWu) 



du 



(Cn.,/i'(i;ja(F.JAW,) 



' du 



du 



lie 



"IIlp(Q;R'') 



\\f^'{Yr,Ja{Yr,JAWu\ 



(3.49) 



which implies that for p > 4 and h E (0, 1] there exists a suitable constant Cp^T such that 

5p — 4 



^42 < 2C;^r/l^ < Cp,T/l'. 

Gathering (I3.48P and (13.501) we derive from (I3.45P that for p > 4 



(3.50) 



(3.51) 
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Similarly, we split J5 as follows: 

k=0 



sup V / (Xt^ - Yk,fi'{Yk)a{Yk)AwAdu 

0<s<t I j-i Jt,. \ ' 



+ 2 


1 f 
sup / 




0<s<t 1 Jt„ 



p 



(3.52) 



With regard to J51, following the same line as f l3.46p yields 



J51 <- 



P-2 ^ JU 



< sup - Ys 

0<s<t 



1/2 



<- sup 

4 o<s<t 



<- sup 11X5 - 

4 0<s<t 



fc=o "'^'^ ' (3.53) 

E ^ / ll/i'(n)a(n)AW-„||i,(^^«,)rf« 



(^2)2 



fc=0 



where the fact Yt^ = Y^, k = 0, 1, X — 1, elementary inequalities and fl3.39p were used. For 
J52, we follow the same way as fl3.49p to obtain 

< sup \\X, - n 



7 \ 2 

o'52 \ 2 ^, £_i 



0<s<t 



2 h~t~^ 



Thus 



J52 < 2 sup ||X,-F,| 

0<s<t 



- 3p-4 1 _ „ 

4 0<s<t ^ ' 



3p-4 



+ 4C;j,/i— . (3.55) 



Note that ^ > 2 for p > 4 and that h<l. Plugging (1333|) and ( ]335|l into ( ]332|l yields 



■h< 7^ sup - + Cp^rh^. 

^ 0<s<t 



■s\\Lp{Q-M'') I 

Inserting ([332]), (1333]), ([331]), f l33T]) and fl336]) into (13:37]) leads to 

2 

LP{Q;IR'*) 

Hence, using estimates in Lemma 13.31 and Lemma 13.61 we derive from (I3.36P that 







/ sup 


Xr — Yj. 


Jo 0<r<u 





au -\ — sup 

LP(n;Rd) 2 o<s<t 



Xs — Yo 



(3.56) 



+ Cp^rh^. (3.57) 



sup \\Xs - Ys\ 

0<s<t 



LP(C;M) 



<2{K + l + mK'' +p'm'K + ^mp'K'') j ^sup \\Xu -Yj(l,^^.^,.^ds 



0<«<s 



+ 7; sup 

^ Q<s<t 



Xs-Y, 
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LP(Q;M'') 



+ Cp^Th''. 



Thus II supo<,<j \\Xs-Ys 



ILp(Q;] 



is finite by Lemma 13.41 and 



sup \\Xs - Ys 

0<s<t 



LP{n;R) 



< a 



sup - Yu\ 

0<u<s 



2 

LP{n;l 



ds + Cp^xh^- 



(3.58) 



The Gronwall inequahty gives the desired resuh for p > 4. Using Holder's inequality gives the 
assertion for 1 < p < 4 and the proof is complete. □ 



4 An illustrative example 

In [S] , the authors have demonstrated the computational efficiency of the tamed Euler scheme, 
compared to the implicit Euler method. In this section we compare computational efficiency of 
the tamed Milstein scheme and the tamed Euler scheme. To this end we choose a simple SDE 

dXt = -Xfdt + XtdWt, Xo = 1 (4.1) 

for t G [0, 1]. Figure [1] depicts the root mean-square errors (11. 3p as a function of the stepsize 
h in log- log plot, where the expectation is approximated by the mean of 5000 independent 
realizations. As expected, the tamed Milstein scheme gives an error that decreases proportional 
to h, whereas the tamed Euler scheme gives errors that decrease proportional to . To show 
the efficiency of the tamed Milstein method clearly, we present in Figure[2]the root mean-square 
errors of both methods as function of the runtime when G {2^°, 2^^} and the mean of 1000 
independent paths are used to approximate the expectation in (II. 3p . Suppose that the strong 
approximation problem (II. 3p of the SDE (14. ip should be solved with the precision e = 0.001. 
From Figure [2l one can detect that = 2^*^ in the case of the tamed Milstein method (ll.lOp 
and that A^ = 2^^ in the case of the tamed Euler method (I1.5P achieves the desired precision 
e = 0.001 in (II. 3p . Moreover, the tamed Milstein scheme requires 8.1860 seconds while the 
tamed Euler scheme requires 147.9230 seconds to achieve the precision e = 0.001 in (II. 3p . The 
tamed Milstein method is for the SDE (14. ip with commutative noise thus much faster than the 
tamed Euler method. 
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